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Introduction: Polar layered deposits are exposed 
in the walls of the troughs cutting the north polar cap of 
Mars. They consist of alternating ice and dust layers or 
layers of an ice-dust mixture with varying proportions 
and are found throughout the cap. Layers thickness 
ranges from meters to several tens of meters with a ~ 
30 meter dominant wavelength [1,2], Although their 
formation processes is not known, they are presumed to 
reflect changes in ice and dust stability over orbital and 
axial variations. Intensive 3-D LMD GCM simulations 
of the martian water cycle have been thus perfomed to 
determine the annual rates of exchange of surface ice 
between the northern cap and tropical areas for a wide 
range of obliquity and orbital parameters values.These 
rates have been employed to reconstruct an history of 
the northern cap and test simple models of dust-ice lay- 
ers formation over the last 10 Ma orbital variations. 

Model description: We use the 3-D water cycle 
model simulated by the 3-D LMD GCM [3,4] with an 
intermediate grid resolution (7.5° longitude x 5.625° 
latitude) and 25 vertical levels. The dust opacity is con- 
stant and set to 0,15. No exchange of ice with regolith 
is allowed. The evolution of the northern cap over 
obliquity and orbital changes (eccentricity, Longitude 
of perihelion) has been recently described with this 
model in [5,6], High summer insolation favors transfer 
of ice from the northern pole to the Tharsis and Olym- 
pus Montes, while at low obliquity, unstable equatorial 
ice is redeposited in high-latitude and polar areas of 
both hemisphere. The disappearance of the equatorial 
ice reservoir leads to a poleward recession of icy high- 
latitude reservoirs, providing an additional source for 
the cap accumulation during each obliquity or orbital 
cycle. Furthering the efforts in [5], a quantitative evo- 
lution of ice reservoirs is here investigated for vari- 
ous astronomical conditions. 

Polar ice annual loss rates (PIALR): the sensivity 
of the annual polar ice rates has been determined from 
12 equilibrated simulations at various obliquities and 
orbital parameters leading to an unstable northern cap 
[5]. Because the most intense cap sublimation occurs 
at the beginning of the summer when the permanent 


cap is exposed, we expect that the summer solstice in- 
solation is, at first order, a good forcing function for 
the {PIALR) variations. However, high summer insola- 
tions correspond to shorter summers (summer solstice 
at perihelion, for example). We found that, at first or- 
der, the summer solstice insolation divided by the mar- 
tian orbital speed at summer solstice is well correlated 
with the PIALR. Fig.l illustrates the marked exponen- 
tial correlation between the PIALR and this new inso- 
lation forcing [7], 
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Figure 1: Correlation between the PIALR and the modified 
summer solstice insolation. The best least-square fit (solid line) 
exhibits a critical insolation close to 300 W.m-2 above which the 
northern cap becomes unstable. 

Polar ice annual accumulation rates (PIAAR): Be- 
low the critical insolation, equatorial ice reservoirs be- 
come unstable. The PIAAR have been determined 
from 6 equilibrated simulations for various obliquities 
(30°, 25.19°, 20° and 15°) and extreme orbital configu- 
rations ( e=0 and e=0.0934 with Lp=90° and Lp=270°). 
Interestingly, the PIAAR have been found nearly un- 
sensitive to obliquity and orbital changes and range be- 
tween 1.4 and 2.0 mm/yr for decreasing obliquities. 
These rates appear also unsentitive to the equatorial 
reservoir size. This can be explained by the fact that, in 
these conditions, most of the equatorial ice is deposited 
first in high-latitude areas. In this context, the PIAAR 
have been considered as a simple function of the obliq- 
uity. 
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Recent Evolution of the Northern polar cap: 

Considering the previous annual accumulation/loss ice 
rates, the thickness of the polar cap has been computed 
over the recent 10 Ma insolation and obliquity changes 
[8] with a 500 yr step. For each orbital condition, the 
accumulation/loss rate is interpolated from the PiAAR 
and PIALR values. The evolution of the surface reser- 
voir has been divided in three budgetary elements 
shown in Fig.2. The evolution of the tropical water 
(30°S-30°N) has been determined from the previous 
simulations and interpolated for an arbitrary orbital 
condition. The exchange between the high-latitude res- 
ervoirs and the tropical budget is determined to ensure 
the global ice mass conservation at each step. 
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Figure 2: Cartoon illustrating the three-box surface ice res 
en’oirs considered. No exchange is allowed between the Other Sur- 
face Ice Reserx’oirs (OSIR) and the northern cap although this ex- 
change is possible when the equatorial reservoir disappears. This 
choice allows to detect the occurrence of these periods. 

Influence of a dust lag deposit: A common pre- 
sumption is that polar water ice contains a small frac- 
tion of dust, which could form a thick lag deposit pro- 
tecting deeper ice deposits from further sublimation at 
higher insolation periods. To incorporate a such 
mechanism, we considered that, when a thickness// of 
ice has sublimed, the PIALR are reduced by a factor f. 

In our standard model, characteristical values //=10 
m and/=10 have been used. However, the total dust 
thickness is not predicted by our model. 

Results: An exemple of the evolution of the surface 
ice reservoirs is shown in Fig.3, starting with a current 
3 km-thickness northern cap. In this case, polar ice rap- 
idly goes to the equator in the high mean-insolation 
period (5-10 Ma ago) and the onset of the polar cap be- 
gins around 4 Ma during the transition towards a lower 
mean-insolation period (0-5 Ma). As predicted in [5], 
high-latitudes deposits could have been formed during 
this transition, when the equatorial reservoir progres- 
sively disappears. We predict that a maximal global 
equivalent layer (GEL) of 100 m could have been de- 
posited in high-latitude areas. After this transition, the 
northern cap grows without major erosional periods 
and each “oscillation” corresponds to the formation of 
a dust-ice layer. 



Figure 3: Evolution of the surface ice reservoirs thickness 
following the three-box cartoon described in Fig.2 and for our 
standard dust lag model. The thickness indicates the equivalent 
ice thicbiess uniformly distributed on the whole surface of the res- 
ervoir. For more clarity, the ice thicknesses of the tropical reser- 
voir and of the OSIR have been multiplied respectively by 100 and 
10 . 

About 620 meters of ice accumulates from the 
equatorial ice reservoir. A statistical analysis shows 
that this corresponds to the formation of 27 ice-layer 
with a 22.0 m averaged thickness and a 13.8 m stan- 
dard deviation. These properties are not changed for 
moderate changes in H and / parameters. This is very 
close to the thickness observed in polar layering [2,8]. 
No accumulation of polar ice from equatorial ice oc- 
curs during the recent low-amplitude obliquity varia- 
tions (0-500 ka and around 2.5 Ma ago), while tran- 
sient equatorial reservoirs are predicted during the re- 
cent summer insolation peaks (0.5-2. 2 Ma). 

Discussion: Our standar model suggests that most 
of the present northern cap thickness does not come 
from the contribution of past equatorial ice reservoirs, 
but more likely from the poleward recession of north- 
ern icy high-latitude deposits during the intervals (0- 
500 ka and around 2.5 Ma). This has to be compared to 
geological observations. However, other surface (South 
polar cap) or subsurface (regolith) ice reservoirs may 
have participated in its formation. 

Critical insolation values are also probably very 
sensitive to surface (polar albedo, thermal inertia) and 
atmospheric properties (atmospheric dust content, C02 
pressure..) and must be investigated in details. This 
could significantly changes the surface ice reservoirs 
history. 
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